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Abstract 

We show the successful intercalation of large area (1 cm^) epitaxial few-layer graphene grown on 
4H-SiC with FeCla. Upon intercalation the resistivity of this system drops from an average value of 
~ 200 Q/sq to ~ 16 Q/sq at room temperature. The magneto-conductance shows a weak localiza¬ 
tion feature with a temperature dependence typical of graphene Dirac fermions demonstrating the 
decoupling into parallel hole gases of each carbon layer composing the FeCla intercalated structure. 
The phase coherence length (~ 1.2/rm at 280 mK) decreases rapidly only for temperatures higher 
than the 2-D magnetic ordering in the intercalant layer while it tends to saturate for temperatures 
lower than the antiferromagnetic ordering between the planes of FeCla molecules providing the first 
evidence for magnetic ordering in the extreme two-dimensional limit of graphene. 
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The chemical functionalization of graphene with atomic species and/or molecules offers 
a unique way to engineer novel physical properties which are not found in pristine graphene 


Q. 


For example the functionalization with fluorine or hydrogen is known to open an energy- 


gap whose value is a function of the coverage of adatoms 


, 3 - 


On the other hand, the 


intercalation of guest atoms or molecules into few layer graphene grants access to unique 
conductive jd], magnetic 3 ^^^1 superconducting joj properties of these materials. Recent 
experiments on intercalation with FeCla of mechanically exfoliated few-layer graphene have 
reported an extraordinarily large charge doping (up to 9 * with a record low 

electrical resistivity (< 8 O/sq) and high optical transparency in the visible wavelength (> 84 
%) which make this material attractive for applications requiring transparent electrodes. 
At the same time the experimental realization of Ca- and FeCls- intercalated bilayer graphene 
a, l6| provides a promising platform to investigate magnetism and superconductivity in two- 
dimensional systems. 

So far most of the scientific effort has focused on the study of chemical functionalization 
of mechanically exfoliated graphene and of graphene grown by chemical vapour deposition, 
whereas the functionalization of epitaxially grown graphene on SiC is still in its infancy 
7|, [Sj . Since epitaxially grown graphene on SiC is one of the main contenders for large area 
manufacturing 3 of high quality graphene, understanding the chemical functionalization 
of epitaxial graphene is a fundamental hrst step to ensure the full exploitation of this ma¬ 
terial. Furthermore, SiC-based electronic devices are uniquely suited for complementary 
applications to those targeted by standard Si-based electronics [l^- Therefore, functional¬ 
ized epitaxial graphene on SiC could enable the development of a new class of electronic 
applications for hostile-environment working conditions -e.g. space applications with high 
radiation levels [h| or nuclear power centrals. From a fundamental point of view, function¬ 
alized epitaxial graphene could enable the study of states of matter such as magnetism and 
intrinsic superconductivity in the extreme two-dimensional limit inherent in this atomically 
thin material. 

Here we present a detailed study of FeCls intercalated few-layer epitaxially grown 
graphene on 4H-SiC with a surface area of ~ 1 cm^, see Figure la. From the analysis 
of the Raman spectra and conductivity measurements we hnd that the intercalation process 
is uniform and the typical room temperature square resistance drops from 174 ±9 VL/sq 
to 16.6 ± 0.6 VL/sq upon intercalation. The measured magneto-conductance (MC) shows a 
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strong weak localization (WL) feature at cryogenic temperatures (< 25K) whose tempera¬ 


ture dependence is well-described by the WL theory of a single layer graphene 


121. These 


hndings demonstrate that the intercalation of FeCla originates parallel single layer-like hole 
gases in the stacking with a phase coherence length (L^) as large as 1.17 ± 0.08/im at 
280 mK. The temperature dependence of shows a steep decrease of for temperatures 
higher than ~ 30iL, a temperature compatible with 2-D magnetic ordering in the FeCla 
layer, while tends to saturate at temperatures lower than ~ 4iF, a temperature compati¬ 
ble with 3-D antiferommagnetic coupling between planes of FeCls. These observations show 
that FeCla intercalated epitaxial graphene is a unique system to study magnetic ordering 
in the extreme two-dimensional limit of graphene. Finally the intercalation of FeCls offers 
a simple way to induce high charge carrier density in graphene (> 1.5 x lO^^cm”^) which is 


pivotal for pioneering the physics of plasmons in unprecedented ranges of wavelengths 


13|. 


Pristine few-layer graphene ( 


terminated face of 4H-SiC 


HIi5 


TG) was grown via thermal decomposition of the silicon 


15j. The layer number of pristine graphene was determined by 
Raman spectroscopy. The intercalation of FLG with FeCls is schematically shown in Figure 
lb. This process was performed in vacuum for a total of 24 hours using the two zone vapour 
transport technique described by Khrapach et al. j4|. Raman spectroscopy and electrical 
transport measurements were employed to characterise the structure and macroscopic sheet 
resistance of the intercalated samples. To ascertain the intrinsic electrical properties of this 
material we measured the magneto-resistance with standard lock-in technique in Hall bar 
devices (length 200 /xm and width 50 /xm) from 200K down to 280mK and with perpendicular 
external magnetic helds up to ± 2T. 

Recent experiments on FeCls intercalated FLG obtained from mechanically exfoliated 
natural graphite have reported a large charge transfer between the intercalant and graphene 
This modihes the Raman spectra of FLG in two distinct ways: an upshift of the G-band 
and a change of the 2D-band from multi- to single-peak structure, respectively (see Figure 
Ic). A graphene sheet which has only one adjacent FeGls layer gives rise to the Gi-band 
(1612cm“^), whereas a graphene sheet sandwiched between two FeGls layers is characterised 


by the G 2 -band (1625cm 


-l^ 


!, 3 [i^ . The changes observed in the shape of the 2D-band 
indicate decoupling of the FLGs into separate monolayers due to the intercalation of FeGls 
between the graphene sheets. 

Figure Id shows a comparison of the measured Raman spectra for pristine and interca- 
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FIG. 1. Panels (a) and (b) show a picture of the intercalated FLG on 4H-SiG and a schematic crystal 
structure of trilayer FeCla-FLG, respectively. The graph in (c) is a plot of the G and 2D Raman 
bands for mechanically exfoliated pristine (top) and FeCls intercalated trilayer graphene (bottom) 
on Si02. d) Shows the Raman spectra for epitaxial pristine FLG (top) and FeCls intercalated FLG 
(bottom). The red curve is a fit of the 2D peak to a multiple Lorentzian function (green curves). 
Upon intercalation with FeGls a single Lorentzian is the best fit of the 2D peak (red). 


lated FLGs. Unlike graphene deposited on Si02, the Raman spectra of epitaxial graphene 
contains the signal from the SiC substrate in the spectral region of the G band. To separate 
the 4H-SiG and graphene Raman spectra we have subtracted the spectrum of 4H-SiG in the 
G-band region by following the procedure introduced by Rohrl et al. 18|. For pristine FLG 
on 4H-SiG, both G and 2D bands are upshifted as compared to pristine FLG on Si02 due 


to the substrate induced strain 


18 


201 1 . To identify the number of layers in pristine FLG 


on SiG we used the method developed by Lee et al, where it was shown that the value of 
the full width at half maximum of the 2D band (FWHM(2D)) exhibits a linear relationship 
with the inverse number of layers (N): FWHM(2D) = (—45(1/N) + 88)[cm“^] [^. In our 
epitaxial FLG samples, the 2D band has a FWHM of 72 cm“^, from which we estimate 
a layer number of three (see Figure Id). Upon intercalation with FeGls, the G-band of 
epitaxial FLG upshifts to the Gi (1610 cm“^) and G 2 (1622 cm“^) positions demonstrating 
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FIG. 2. a) Shows a plot of the measured four terminal resistance (R) as a function of the ratio 
between the conductive channel lenght (L) and width (W). The blue data points correspond to 
pristine FLG on 4H-SiG while the red data points correspond to FeCla-FLG on 4H-SiG. The inset 
shows a scheme of the configuration used to measure the resistivity, b) and c) are histograms 
showing the distribution of resistivity for pristine FLG (Blue) and FeCls-FLG (Red), respectively. 

the successful intercalation of FeCla molecules in the epitaxial FLG. At the same time, the 
2D-band structure of the FLG on SiG is changed from a multi-peak structure to that of 
a single Lorentzian and the 2D-band downshifts to that of monolayer graphene on 4H-SiG 
(«i2715cm“^), see Figure Id. These observations demonstrate that our FeGla-FLG on SiG 
has three graphene sheets and two FeGla layers intercalated between the graphene sheets 
as schematically illustrated in Figure lb. Therefore, this intercalation stage leads to the 
formation of 3 parallel hole gases, where the top and bottom gases will have similar values 
of charge density which will be lower than the charge density of the middle hole gas. 

The sheet resistance of the as grown epitaxial FLG and FeGls-intercalated epitaxial FLG 
was hrst measured in a four terminal conhguration in macroscopic samples with a conductive 
channel of hxed width (0.7 cm) and channel length ranging from 0.7 mm to 4.2 mm using 
Ti/Au (5/50 nm) contacts. It was found experimentally that the pristine epitaxial FLG 
have an average sheet resistance {Rsq) of 174 ± 9 G/sq, estimated from the linear £t shown 
in Figure 2a. The corresponding values of resistance in the same devices after intercalation 
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with FeCls are systematically much lower with an average value of 16.6 ± 0.6 hl/sq. The 
distribution of Rgq measured for FLG and FeCla-FLG are presented as two histograms, see 
Figure 2b and c. In both cases we find a single peaked distribution with narrow spreading 
of the values of Rgq demonstrating the homogeneous and reproducible electrical properties 
in large area FeGla-FLG. 

The charge carriers sign and density are readily characterized from measurements of the 
transverse magneto-resistivity (pxy) in Hall bar devices (see inset of Figure 3a). The positive 
slope of pxy(B) stems for holes as the dominant charge carrier type with a typical density 
of > 1.5 X 10^^ cm“^. This high charge carrier density far exceeds the values achievable 
with top-gates on SiG|2]| and is comparable to that achieved in electric-double-layer FLG 


transistors 


22 


23|. 


Measurements of the temperature dependence of the longitudinal resistivity reveal that 
Pxx decreases monotonously from 175K down to 25K as expected for heavily doped graphene 
J], see Figure 3b. However, at temperatures below 25K we systematically observe an in¬ 
crease in pxx which was not previously observed in intercalated mechanically exfoliated FLG 
and whose quantum nature can only be unveiled with a detailed study of the magnetic field 
dependence. Figure 3c shows a typical measurement of Pxx(B) at teinperature T=300mK. 
A zero-held peak characteristic of weak localization is clearly visible 2^. As the magnetic 
held breaks time-reversal symmetry, the quantum interference is reduced and a classical 
magneto-resistive signal is observed (|B| > ±250mT). Note that owing to the macroscopic 
size of the epitaxial-FLG Hall bars we can directly measure the weak localization signal 
without the need for subtracting any mesoscopic conductance huctuations typically present 
in devices fabricated on mechanical exfoliated graphene 241. 

A detailed study of the temperature dependence of the longitudinal magneto-resistivity 
shows that the WL peak is heavily suppressed when the temperature increases, see Figure 
4a. At the same time the concavity of the roughly parabolic shape background of p^x vs. B 
is clearly temperature independent. Recent experiments have shown that the longitudinal 
magneto-resistance can be understood as the sum of a Drude term plus a correction due to 
electron-electron interaction. The hngerprint of the electron-electron interaction correction is 
a strong temperature dependent concavity of Pxx(B) [25|. The fact that our measurements 
show no signihcant change of the concavity as a function of the temperature (see Figure 
4a) implies that for FeGls-FLG the electron-electron interaction correction to the magneto- 
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FIG. 3. a) Shows a plot of the Hall resistivity {pxy) measured in a representative Hall bar device 
shown in the micrograph picture in the inset, (b) and (c) show the temperature and perpendicular 
magnetic field dependence of the longitudinal resistivity {pxx) measured in the same device. 


resistance is negligible compared to the classical Drnde term. Indeed, if we only consider 
the condnctivity tensor with terms (Jxxj = and cXxyj = where aoj is the 

classical Drnde condnctivity for each hole gas indexed by j, we can describe accnrately the 
experimental data nsing the eqnivalent magneto-resistance for 3 parallel hole gases of which 
two with similar charge density as previonsly identihed from Raman spectroscopy, see Fignre 
4a. 


Farther evidence that the hole gases present in FeCla-FLG are deconpled Dirac fermions 
is provided by a detailed analysis of the WL correction which we condnct on the signal 
after snbtracting the classical magneto-resistance backgronnd, see Fignre 4b. We consider 
the eqnivalent MC for the intercalated-FLG to be Aa = Aa^, with Acr^ the well- 


established MG for Dirac fermions in single-layer graphene 


121 : 


Aa,{B) 



-F 


B 




-2F 


B 


B^j + Bij + 


( 1 ) 


7 











B is the applied perpendicular magnetic field, and diffusion coefficient 

D = hvF^/n/2e^pxx'/n, F{z) = Inz + {\ + \') and tjj is the digamma function. 

In graphene WL is uniquely sensitive to inelastic phase breaking processes characterized 
by a dephasing scattering time and also to elastic scattering that breaks the chirality of the 
charge carriers (e.g. scattering off surface ripples, dislocations and atomically sharp defects) 
with a corresponding scattering time A potential anisotropy of the Fermi surface in k 
space, i.e. trigonal warping, can also destroy intravalley WL and the associated scattering 
time is and Finally the intervalley scattering which occurs at a rate 

on defects with size of the order of the lattice spacing, promotes the interference of carriers 
from different valleys in momentum space. 

Figure 4b shows the measured WL signal extracted from the longitudinal magneto¬ 
resistivity after subtracting the classical background and the corresponding best fit of the 
total MC modelled with the parallel of three decoupled hole gases for different tempera¬ 
tures. At T=280mK the best-fit values for the middle hole gas with higher charge density 
are Bi = 7.4 * 10“^ ± 6 * 10“^ T and B^ = 1.2 * 10“"^ ± 2 * 10“^ T and for the outer hole 
gases with lower charge density Bi = 4.6 * 10“^ ± 5 * 10“^ T and B^ = 1.6 * 10“^ ± 3 * 10“^. 
In this fit procedure we find that the term in Eq. (1) containing B^ is negligible, stemming 
for a very short r* compared to the other characteristic times (i.e. r* and r^). The good 
agreement between the experimental data and the model holds for a very wide range of 
temperatures starting from 280mK up to the highest measured (80K). We can therefore 
conclude that our samples are in the regime where Ttr ^ Ti, « Ti « t^. If we assume an 
identical diffusion constant for each hole gas we can estimate the value of at each different 
measured temperature. Since the typical behaviour of the dephasing time, as predicted by 
theory is power-like 261] (e.g. oc T~‘^F for narrow wires), we plot our data in a double 
logharitmic scale, see Figure 4c. This analysis reveals two different regimes: rapid decrease 
of the dephasing time for temperatures T > AK and nearly saturated behaviour at low tem¬ 
perature. Such saturation is common for samples with magnetic impurities in microscopic 
concentrations j^, ^ which in graphene can originate from chemisorbed magnetic atoms, 
hydroxyl groups, vacancies or resonant localized states in the vicinity of the Dirac point 29]. 
In our case the presence of iron in the intercalant specie is a source of magnetic impurities, 
which can account for the observed saturation of the dephasing time. 

The dephasing length can be directly extracted from the parameter B^ using the relation 
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FIG. 4. The graph in (a) is a plot of the measured magneto resistance (scatter points) at various 
temperatures (see colour coded legend). The continuous lines are a fit to the Drude conductivity, b) 
Shows the measured weak localization contribution to the magneto-conductance (scatter points) 
and the best fit (solid line) to the model of Eq. (1) for various temperatures, c) shows the 
temperature dependence for the extracted for two different intercalated samples (labelled SI 
and S2) and (d) shows a summary of the temperature dependence of the values of present in 
literature relative to pristine single-layer graphene grown and/or obtained by different methods as 
specificed in the legend and a comparison to the values of estimated in our FeCla-FLG. 


4eLj 


Therefore without making any assumption on the 


which gives = 

diffusion constant of the hole gases we extract the corresponding values of for a wide range 
of temperature (from 80K down to 280mK), see Figure 4d. Upon lowering the temperature 
the values of estimated for FeCla-epitaxial FLG increase monotonically, however below 
~ AK this increase is less pronounced with a maximum value of = 1.17 ± 0.08/rm 
at 280mK. To understand the possible physical origin of the temperature dependence of 
Lfj, in our material, we compare L^{T) for two representative FeCLa-epitaxial FLG (i.e. 


indicated by SI and S2) with values measured in pristine epitaxial FLG 


25 
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GVD 
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graphene 3l| and graphene mechanically exfoliated from natural graphite 32|]. The overall 
temperature dependence of is similar for all the different kinds of graphene, however 
two striking features distinguish FeCla-epitaxial FLG from pristine-FLG. More specihcally 
FeGls-FLG has (1) larger values of than pristine graphene at a given temperature and (2) 
Lfj, decreases rapidly only for temperatures much higher than what is observed in pristine 
graphene, see Figure 4d. 

To explain the macroscopic scales of the quantum coherence estimated in our samples even 
in the presence of magnetic FeGls molecules - expected to be a source of excess dephasing- 
and the temperature dependence of we have to consider what is known on the magnetic 
ordering in stage one FeGls-intercalated graphite. This material is a Heisenberg system 


characterized by two distinct magnetic orderings 


33 


35| . these are (1) strong 2-D magnetic 


correlations in the plane of FeGla with a transition temperature of ~ 30K and (2) a tran¬ 
sition to 3-D antiferromagnetism between planes of FeGb at 3.8K. These magnetic phases 
lave been studied experimentally in bulk materials with Mossbauer and neutron diffraction 


33l. l34| . A theoretical understanding of the structural, electronic and magnetic properties 


in FeGls-FLG was only recently presented using ab initio numerical simulations using a 
generalized gradient approximation of Perdew-Burke-Ernzerhof with an effective interaction 
oarameter which accounts for the Hubbard parameter as well as the exchange interaction 
361]. It was found that for stage one FeGls-FLG within the spin-polarized calculations at 
low temperature the 3-D magnetic ground state is antiferromagnetic, that is 1.5 eV per unit 
cell lower than the ferromagnetic conhguration. At the same time the magnetic moments 
of FeGls are expected to align ferromagnetically in the same intercalant layer. Finally, the 
band structure of this compound is also expected theoretically to have a linear dispersion 
with a Fermi energy of ~ 0.9eV 16, 36] as conhrmed by recent experimental studies jl]. 


Spin-flip scattering off time dependent ranc 


of excess dephasing in mesoscopic systems 27H29||. However, upon attaining a macroscopic 


omly oriented magnetic moments is a cause 


magnetic ordering, the charge carriers experience less spin-flip scattering events leading to 
a longer phase coherence length. Although the material studied in this article is not bulk, 
i.e. it is only a 3-layer FLG, we observe experimentally large values of L^p for T < 30K 
(compatible with the 2-D magnetic correlations in the plane of FeGls) and a saturation of 
for T lower than ~ 4K, see Figure 4d. Therefore, considering what is known on the 


magnetic ordering in this kind of intercalated compounds 


33 


36| . our experimental Endings 
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strongly suggest that magnetic ordering of the FeCls planes suppress magnetic scattering in 
FeCls-FLG. For temperatures higher than the 30K a sharper decrease of is observed in 
FeCla-FLG as compared to pristine graphene indicating that randomly oriented magnetic 
moments in the intercalated FLG are driving excessive dephasing. 

In conclusion we have demonstrated the successful intercalation of large area epitaxial 
graphene grown on 4H-SiG with FeGla. We found that this material has a very low square 
resistance (~ IGhl/sq) becoming an ideal candidate for transparent electrodes in novel com¬ 
plementary applications targeted by SiG-based electronics. With a detailed study of the 
Raman spectra we have characterized the structure of this material. Finally we have con¬ 
ducted a systematic study of MG in Hall-bar devices and we have estimated a large dephasing 
length of Lrj) ~ 1.2/im at 280mK in FeGls-FLG. The observed temperature dependence of 

is compatible with the occurrence of magnetic ordering in stage one FeGls-FLG, showing 
that this novel material is a good platform for studying magnetic ordering in the extreme 
two-dimensional limit of graphene. 
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